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Relation between the infrared renormalons, the Borel resummation prescriptions, and the analyt¬ 
icity structure of Green functions in perturbative QCD (pQCD) is investigated. A specific recently 
suggested Borel resummation prescription resulted in the Principal Value and an additional power- 
suppressed correction that is consistent with the Operator Product Expansion. Arguments requiring 
the finiteness of the result for any power coefficient of the leading infrared renormalon, and the con¬ 
sistency in the case of the absence of that renormalon, require that this prescription be modified. 
The apparently most natural modification leads to the result represented by the Principal Value. 
The analytic structure of the amplitude in the complex coupling plane, obtained in this way, is 
consistent with that obtained in the literature by other methods. 
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Green functions in QCD are quantities which, in general, are known to possess renormalons, i.e., singularities of 
the Borel transform on the real axis [jj. They appear as a consequence of the specific asymptotic behavior of the 
high-order perturbation coefficients. The singularities on the positive axis, called infrared (IR) renormalons, represent 
an obstacle to the Borel integration and lead to the well-known renormalon-induced ambiguity for the observable. 
The most widely used prescription for fixing this ambiguity has been to perform the Borel integration parallell to 
the positive real axis and taking the real part of the resulting integral (Principal Value; see, e.g., Refs. ptyj). This 
prescription in pQCD was shown to be favored by or to be consistent with analyticity requirements in the momentum 
plane |4| and in the coupling parameter plane |5J. Further, this prescription is attractive due to its mathematical 
simplicity. If the QCD quantity under consideration allows to be presented by the Operator Product Expansion 
(OPE), then any additional, nonperturbative, terms are expected to have the power-suppressed form of the higher- 
twist terms. The perturbative QCD (pQCD), involving the quark and gluon degrees of freedom, is usually expected 
to be unable to predict the strength of such terms. 1 

Recently, a specific prescription, based on the IR renormalon considerations, has been proposed @ to fix the strength 
of such higher-twist terms, and this method gave encouraging numerical results in the case of the resummation of 
the Gross-Lewellyn-Smith sum rule and of the heavy-quark potential Jt] . The main observation was that the IR 
renormalon induces in the Borel-integrated quantity a nonphysical cut along the positive axis in the complex plane 
of the coupling parameter z, and that this cut structure can be naturally eliminated by subtracting a cut function 
proportional to (—z)" where v is related with the power coefficient of the renormalon singularity. Further, the energy 
dependence of the subtraction term, when the coupling parameter z is positive, is consistent with the predictions of the 
OPE for the corresponding higher-twist term. For positive z, the result is the Principal Value of the Borel integration 
minus the aforementioned term. A somewhat speculative interpretation of this result suggests that in this way the 
dominant part of the genuine nonperturbative higher-twist effect is obtained (the correction term to the Principal 
Value), although the method is based on perturbative (pQCD + renormalons) knowledge only. A more conservative 
rephrasing of this would be that this result represents “the most that we can get” out of pQCD, i.e., the natural 
basis to which one should eventually add other contributions to the aforementioned higher-twist term; such genuine 
nonperturbative contributions would involve the vacuum expectation values of the higher-twist operators appearing 
in the OPE. 

In the present work, the aforementioned method is scrutinized. As a result, more support is given to the second of 
the two mentioned interpretations of the pQCD renormalon resummation results. Even more so, the results suggest 
that the natural resummed perturbative contribution is just the Principal Value of the Borel integral, without the 
aforementioned higher-twist term. As a by-product, some insights into the role of the IR renormalons in the analyticity 
structure of the pQCD amplitudes in the coupling plane are obtained. 

Let us consider a Euclidean QCD amplitude A[a(Q)] where the quark mass effects are neglected. Therefore, it 
can be regarded as depending on the energy Q = \/—q 2 of the corresponding process only via the QCD coupling 
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parameter a(Q) = a s (Q ; MS)/ 7 r whose running is determined by the (MS) renormalization group equation 


ttt 4 = Po a 2 (1 + C! a + c 2 a 2 + c 3 a 3 -\ -) . (1) 

am /j, z 

For simplicity of argument, we will consider only the effects of the leading IR renormalon and will neglect the 
subleading IR renormalons. We will also assume that such an amplitude can be described by the OPE. 

First, let us review the method of Ref. @ in detail. The presentation of the method is here somewhat different. 
The Borel transform B(b) of A[a(Q)] has a singularity at b > n (n = 1, or 2, or 3, ...) of the form 


B{b) 


C 

(1 - b/n) 1+u 


[l + K\ (1 — b/n) + K 2 (1 — b/n ) 2 • • • ] + analytic part . 


( 2 ) 


Here v = (nc\— 72 n)//? 0 j where 72 n is the one-loop coefficient of the anomalous dimension of the corresponding higher- 
twist operator (with dimension d = 2n) in the OPE for A[a(Q)]. In the following, we will ignore the contributions 
of the terms with coefficients Kj (j > 1 ), since their inclusion is straightforward (zz i—> v — j) and does not affect the 
conclusions. For definiteness, the renormalization scale /r is taken to be fi = Q, and the renormalization scheme is 
also regarded as fixed, e.g., MS. We will regard the coupling to be complex in general, z = /3oa(Q)/n= \z\ exp(*</>), 
corresponding to complex momenta Q. Then the Borel integral Abi(z) can be defined as 


Ab,<2) = A 


+ 00 exp(i 0 ) / ^ 

db exp I-) B(b) , 


nz 


(z = |z|e^, 6 = | 6 |e i?i ) . 


( 3 ) 


For z near the positive real axis z = \z\ ±ie [\z\=p 0 a(Q)/n > 0], it is straightforward to show from 0 , with the help 
of the Cauchy theorem, the following formula: 


Abi(*)L = 


z=\z\±ie 


1 


r +oo±ie / b 

db exp (-— I B{b) (s —> +0) . 


'±4£ 


n\z\ 


( 4 ) 


The real part of this is the Principal Value. As a consequence of the IR singularity 0 at b > n, the integral 0 has 
a discontinuity (cut) at the positive real axis z > 0. Namely, the quantity 0 is not real and its discontinuity shows 
up in its imaginary part. This can be seen by introducing in the integrand of Eq. 0 the new complex integration 
variable t via b = n( 1 + |z|f) 


ImAeifj - Wife) = ^ l dbexp (- A.) ^ ^ P (-0 si„(™) W - . (5) 

The (Hankel) contour C t is depicted in Fig.HJ and the expression 0 is obtained from the known Hankel contour form 
of the Gamma function (see, for example, | 8 |) 

T(s) sin(Trs) = ^ = +^t J c dte"*(-t)- 1+s (|s| < oo) , ( 6 ) 


in the special cases s = —v,—v+ 1,.... Because (— |zj =F*e) = \z\ exp(=Fi 7 r) and thus (— |z| =F*e) " 
*sin( 7 ri/)], Eq. © leads to 

Abi(^) = er x ! z r(-i/) (-z) _l/ + A B i(z), 

Po 


|z| ^[ 008 (^ 2 /) ± 

( 7 ) 


where Abi(z) is a function without cuts in the complex z-plane since the first expression on the right-hand side 
absorbs the renormalon-induced cut © at z > 0. The first expression, when z= |z|±ie = /?o a(Q)/n±ie is at the real 
positive axis, has the Q-dependence Q~ 2n a s {Q)' l2n -^ 0 [\ + O{a s )], the same as the corresponding power-suppressed 
(higher-twist) term (02n)^/ (Q 2 ) n in the OPE. The imaginary part of this expression, i.e., expression 0, must be 
identified as the imaginary part of the contribution from the leading IR renormalon ©. The central assumption of 
the method is that the full first expression on the right-hand side of 0 , which contains the full nonphysical cut © 
along z > 0 and no cut along z < 0 , represents the nonphysical cut-function which is to be eliminated 


A (cuti)(*) = -C— e 


~ x ' z F(—i/) ( -z)~ v = C^- e " 1 / 2 
Po 


r(l+i/) sin(7r^) 


(-*)- 


( 8 ) 


This same cut contribution, for z = |z|±*e at the positive real axis, can be obtained also by Borel-integrating the 
nonanalytic part of the Borel transform 0 along its cut b > n above or below the real axis in analogy with the full 








FIG. 1: The path Ct in the integrals J^J and Ihl . 


Borel-integrated quantity 0 . To show this, we use the new real integration variable t such that b = n(l + zt), and 
therefore b > n±ie corresponds to t > 0 

i poo ±ie / 7 \ 

A < g° nan -\z = \z\±ie; b >n) = — / db exp (-(9) 

PO Jn±ie \ nz J 

n f + °° 

= +C— z e~ 1/z / dt e~* {-z)- x - v t~ l ~ v ( 10 ) 

Po J 0 

= -c^ e- x ' z r(-iz) (-*)"*'. (11) 

For z away from the real positive axis, the analytic continuation of this expression in z keeps its form o unchanged, 
i.e., precisely the cut-function ©• The Borel-integration over t > 0 ( 44 > b > n ± ie) in Eq. m converges only when 
Re(i/) < 0 and gives the result (Hill . When Re(z/) > 0 , the result (Hill represents the analytic continuation in v. 

The subtraction of the cut-function contribution 0 leads to the final result for the resummed value of the observable 
at the positive real value z = \z\ = Poa(Q)/n > 0 


A [z = p 0 a(Q)/n] = A B i(|z| ± ie) - A (cutl) (|z| ± ie) = A B i[z = fl 0 a(Q)/n] 

R(b) 


foatie db 

r b 1 

Lie A eXP 

Poa{Q)_ 


(1 - &/n) 1+jy ’ 


( 12 ) 


In the Borel integration here, the exactly known IR renormalon singularity has been factored out explicitly. Function 
R(b) = (1 — b/n) 1+v B{b) (whose truncated perturbation series up to ~ b 2 is known exactly in the case of several 
QCD observables) has a much weaker singularity at b = n than B(b). In Eq. H 2 ll . the equality ReA( cutl )(|z| ± ie) = 
±cot( 7 r^) ImA( cutl )(|.z| ± ie) was taken into account, a direct consequence of Eqs. 0, 0, and (— \z\ =F ie)~ l/ = 
|z| _ 1 "[cos( 7 ri/)±*sin( 7 ri/)]. The results m and m are those arrived at in Ref. . 

The renormalon power coefficient v = (nci — 72 n )/Po should be regarded as a general parameter which can take on, 
in principle, any (real) value. For example, when varying the number of effectively active quark flavors n/ continuously, 
v changes continuously. Yet another example is given by the large-/?o approximation, when ci —> 0 and v — > 0. The 
residue C of the renormalon in general does not vanish in the large -Po limit. 

The cut-function 0 , as a function of general v and when the coupling parameter is near the positive real axis 
z± = \z\ ± ie = Poa(Q)/n ± ie, can be rewritten in the following form: 


A(cuti)(- = \z\ ±ie) 


n_ _i/| z | [~ 7 rcot (7 Tv)Tin\ 

Po 6 r(l+i/) 


(13) 


The central assumption of the method, i.e., the subtraction of the cut-function 0 from the Borel-integrated value, 
appears to be plausible and natural, especially because the cut function © [o m] has a simple form, and because 
it represents precisely the contribution of the Borel integration of the nonanalytic part of the Borel transform 0 
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along the cut b > n parallel to the real axis, as explained in Eqs. © ClUi. However, there are at least two problems 
with this method when we impose on it the plausible condition that it should work for any (real) value of the power 
coefficient v. 

1. When v is nonnegative integer {v = k\ k = 0,1,2,...), the cut-function Cl is infinite, because cot( 7 r^) diverges 
there. This would not present a problem if the residue C of the renormalon disappeared, making this cut- 
function finite. This appears to be unlikely, as argued above; in particular, for Adler function in the large-/?o 
approximation (n = 2, ci = 0 ,74 = 0) this would mean the disappearance of the leading IR renormalon. Note 
that, in contrast to cot( 7 r^), the factor 1/T(1 + v) in Eq. (II dll is an analytic function in the entire complex 
I'-plane. 

2. When v is a negative integer v = —1, —2,...,..., the IR renormalon singularity 0 with the cut disappears (even 
when C 1 / 0), the Borel transform is analytic. This means that A( cutl )(;r) must be zero. However, according to 
Eq. (THU) A {cutl )(z) 7 ^ 0 - note: 7r cot(7r^)/T(l + v) = (—1 ) k {k — 1)! ^ 0 when v = — k = —1, —2,.... Again, it 
is the poles of cot(niy) that cause the problems, but this time at v = —1, —2,.... 

If we take the view-point that the method should be taken as the starting point nonetheless, because of the 
mentioned plausibility and naturalness of the choice of 0 , we should definitely modify it so that the aforementioned 
two problematic points are eliminated. This can be done in a natural way, by inspecting again the expression 
O- The two problematic aspects arose because of the poles of the term cot( 7 r^), at v = 0, ±1,±2, etc. The 
function f{y) = cot(- 7 w) should thus be regularized in order to eliminate both problems. The apparently most natural 
regularization is obtained by subtracting the simple pole functions with the coefficients being equal to the residues of 
f{v) at those poles. This gives: 


— 7TCOt(7w) 


— 7TCOt(7w) -|-h ^ 


fc =1 


v — k v + k 

1 \ 

■ C0t(7 tv) H-h 2v —-— = 0 

7/ • ^ 7/ z — 


1 


(14) 

(15) 


We see that the function cot(7r v) is identical to the sum of the corresponding simple pole functions. Only some of the 
meromorphic functions have this remarkable property. The procedure m d therefore eliminates the real part of 
the cut-function o at the real axis z = |,sj±z£, the modified expression there is purely imaginary. The final result 
of the modified method is then just the Principal Value (PV) of the Borel integral, when a(Q) > 0 


A [z = p 0 a(Q)/n ] 


m 


/3 0 «(O)J (1 — b/n) 1+v 


(16) 


and the new cut-function can be written for a general complex z as 


Acut^ ( z ) = +Cj q e-V» 


r(l + i/) sin(7w) 


[(— z) v — cos( , 7rr / ) z "] . 


(17) 


It can be checked explicitly that this function is finite for any finite complex u, including v = 0, ±1, ±2,.... It is an 
analytic function of z outside the real axis, with the cut along the real axis. The subtraction of the poles of cot(7r^) 
introduced in terms proportional to z ". The consistency of the results for v = 0. ±1,... thus introduced 
an additional cut along the negative real axis. 

In general, any regularization in v of the singular cot(7w) factor in Eq. IIRil) . not just the simplest regularization 
would give an acceptable result. The factor l/T(l+r/) is already nonsingular (even analytic) for all v. The factor 
exp(—l/I^DIzl”" in Eq. ® cannot be modified because it reflects the correct (^-dependence of the aforementioned 
higher-twist term (C? 2 n) (<3 ^/ {Q 2 ) n in the OPE. The imaginary part in Eq. Cl cannot be modified because it is 
needed to make the amplitude A (z) = Abi(-z) — A cut (^) real for the real positive z. Therefore, the most general 
I'-regularization of the cut-function d for z = |z|±*e is represented by a simple substitution cot(m/) 1 —> g(v), where 
g{y) is an arbitrary nonsingular (possibly analytic) function of v which is real for real v 


A cut (.-=M±fe)=-cf 

po r(i+i/) 

This cut function has then the following form in the complex coupling plane: 


■<*> ~ e ~ V * r ( i+»-) 


{(— 2 ) v — [cos(7rv) — g(v) sin(7r^)] z ’'} . 


(18) 


(19) 












5 


Again, as in Eq. G3, we see that the ^-regularization introduces an additional cut along the negative axis z < 0. 
The amplitude A = Abi — A cut is now represented by the Principal Value and a higher-twist term proportional to 
g(y) where g(v) is nonsingular in v. We stress that any such choice is physically acceptable, the Principal Value 
choice being distinguished in this context only by the mathematical simplicity of the corresponding ^-regularization 

o o 

The main idea behind the result Eqs. ©, im as stressed in Ref. was that the cut-function which is to be 
subtracted from a Borel-resummed QCD amplitude with an IR renormalon has a cut only along the positive axis in 
the coupling plane z = a(Q). In the two exactly solvable non-QCD examples presented in Ref. fjj this idea was shown 
to hold. Here we showed that the requirement of finiteness with respect to the IR renormalon power parameter v 
in QCD amplitudes, and their consistency in the absence of the IR renormalon, imply that the cut function must 
contain, in addition to the cut along the positive z-axis, also a cut along the negative z-axis. This, in turn, implies that 
the amplitude A(z) = Abi(z) — A cut (z), while having no cut on the z > 0 axis in accordance with the unitarity and 
causality conditions |fj, does have a cut along the z < 0 axis (Landau regionlas a consequence of the IR renormalon. 2 
This somewhat counterintuitive conclusion was also obtained in Refs, jj, ;,5| by substantially different approaches, 
where the Principal Value prescription was adopted. If the considered amplitude has ultraviolet (UV) renormalons, 
it should have a cut along z < 0 even if it has no IR renormalons. If the coupling a(Q) is regularized so that it is 
finite for all Q 2 > 0 (see, e.g., 0,0, EH 0, Gil), be., in contrast to the pQCD a(Q) it has no Landau singularities, 
the conclusions about the analyticity of the considered amplitude in the complex coupling plane probably change 
singnificantly. 

The IR renormalons play an important role also in the resummations using modified Borel transforms where the 
entire integrand in the Borel integration is renormalization scale (RS) invariant. Such transforms were introduced 
by Grunberg [l7| on the basis of a larger class of transforms proposed in Ref, fl8l in a somewhat different context. 
Such RS-invariant Borel transform resummations were applied in Refs. 0 l20l ]. by either evaluating the Principal 
Value of the Borel integral 0] or adding to the Principal Value the higher-twist OPE terms (20]. The discussed 
method of Ref. for the ordinary Borel transforms 0 , can be adapted to the method of the RS-invariant Borel 
transforms. The problems (divergences) appearing in this case are similar to those discussed here, but algebraically 
more complicated. It is not clear whether in such case an analogous regularization procedure as the one presented 
here would lead naturally to the Principal Value of the RS-invariant Borel resunnnation. 

In the present work, the method of subtracting a power-suppressed term from the Principal Value of the Borel 
integral for QCD amplitudes with IR renormalon, recently proposed in Ref. b] and applied in Refs. ]7], was scrutinized. 
It was pointed out that the result becomes physically untenable for specific values of the renormalon power coefficient 
v, as a consequence of the divergences of a term [cot(7ri/)] appearing in the result. When these divergences are removed 
in apparently the most natural way, the power-suppressed term of the method disappears and the modified result 
becomes the Principal Value. Any removal of the aforementioned divergences results in an IR-renormalon-induced 
cut along the negative axis (Landau region) in the coupling plane. These conclusions suggest, among other things, 
that the most natural pQCD Borel integration of a QCD observable remains the Principal Value. The additional 
power-suppressed (higher-twist, higher-dimensional) terms cannot be inferred from pQCD (+renormalon) methods 
in any natural way. Such additional (OPE) terms involve vacuum expectation values of higher-twist operators and 
can theoretically be obtained or estimated only by genuinely nonperturbative methods. Phenomenologically, such 
additional OPE terms can be determined by fitting them to the corresponding experimental data. However, in such 
a procedure, it is important to keep for the leading-twist term in the OPE a specific resummed pQCD expression, 
most naturally the Principal Value of the Borel integral, i.e., Eq. m- On the other hand, if the leading-twist term is 
taken to be a truncated perturbation series (TPS ) , the strength of the higher-twist terms will sometimes dramatically 
change when the order of the TPS is changed |2ll |22|. 


2 According to Ref. 191. th e physical singularities in confined theories are generated by the physical hadron states. The point 2 : = 0 is an 
essential singularity lldllll . The confinement is not seen by pQCD, thus the cut along 2 : < 0 (Landau region) is not physical but must 
appear in quantities involving pQCD coupling a(Q). 
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